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ABSTRACT 



Aims. Our aim is to assess the benefits and limitations of using the redundant visibility information in regular phased array systems 
for improving the calibration. 

Methods. Regular arrays offer the possibility to use redundant visibility information to constrain the calibration of the array inde- 
pendent of a sky model and a beam models of the station elements. It requires a regular arrangement in the configuration of array 
elements and identical beam patterns. 

Results. We revised a calibration method for phased array stations using the redundant visibility information in the system and ap- 
plied it successfully to a LOFAR station. The performance and limitations of the method were demonstrated by comparing its use on 
real and simulated data. The main limitation is the mutual coupling between the station elements, which leads to non-identical beams 
and stronger baseline dependent noise. Comparing the variance of the estimated complex gains with the Cramer-Rao Bound (CRB) 
indicates that redundancy is a stable and optimum method for calibrating the complex gains of the system. 

Conclusions. Our study shows that the use of the redundant visibility does improve the quality of the calibration in phased array 
systems. In addition it provides a powerful tool for system diagnostics. Our results demonstrate that designing redundancy in both the 
station layout and the array configuration of future aperture arrays is strongly recommended. In particular in the case of the Square 
Kilometre Array with its dynamic range requirement which surpasses any existing array by an order of magnitude. 

Key words, radio astronomy, calibration, redundancy, phased array, phased array calibration, high dynamic range imaging. 



1. Introduction 

An important conceptual difference between the next generation 
of radio telescopes and conventional ones is their hierarchical 
system architecture . An excellent exa mple is the Low Frequency 
ARray ( LOFAR). Ide Vos et al.l(l2009l) . In LOFAR, a station con- 
sists of phased arrays, either sparse or dense. In the sparse 
phased array stations, the station elements are dipoles which are 
digitally beamformed to synthesize a station as a dish such as 
the LBAs (Low Band Antennas) in a LOFAR station. Phased ar- 
rays operating above ~ 100 MHz are often implemented as com- 
pound elements or tiles such as the HBAs (High Band Antennas) 
at the LOFAR stations or in E MBRACE (Electronic Multi-Beam 
Radio Astronomy ConcEpt), lArdenne et al.l (120041) ; iKant et at] 
(1201 lh . EMBRACE is an example of a very dense phased array 
station. In these stations, the station elements are phased array 
tiles (e.g. Fig.[TJ. A tile is a regular arrangement of many dipoles 
whose signals are added in phase to form an instantaneous beam 
(the beamforming). The tile output signals are digitally phased 
up to synthesize a station as a dish. At the next level in the beam- 
forming hierarchy, the beam formed output of each station is 
transported to the central correlator to synthesize the whole tele- 
scope. Calibration has to be performed at different levels of this 
hierarchy to provide a final, high dynamic r ange image of part 
of sky as explained by Wiin holds et al.1 d2010h . 

In this paper, we are concentrate on the calibration at station 
level. Its purpose is to ensure the station beam stability over time 
and frequency. A robust calibration as part of the beamforming 
process should guarantee a stable beam pattern of the station for 




Fig. 1: 24-tile HBA station. Each tile is one station element. This 
is the station configuration for most of LOFAR' s HBA stations 
including the LOFAR core or CS302. It is clear that a station like 
this is highly redundant. 



data going to the central correlator. This is crucial for the dy- 
namic range of the final images made using data from the entire 
array. 

In a phased array station, the output of all station ele- 
ments can be correlated. These correlations are called station 
visibilities which are used for engineering purposes e.g. sta- 
tion calibration, RFI detection/mitigation. These correlations 
include many short baselines on which extended structures 
such as the galactic plane are captured. The most commonly 
used calibration methods for phased array stations are model 
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based, such as a multi-source cali bration method introduced by 
IWiinholds and van der Veenl(l2009al) . A model based calibration 
method requires the presence of one or more relatively unre- 
solved point sources such as CasA and a model of the extended 
structures (see Figj2j. Modeling an extend ed structure is com- 
putation ally difficult and expensive. IWiinholds and van der Veer] 
(2009b) phenomenologic ally model it as correlated noise and es- 
timate the parameters of interest for calibration using a WALS 
(Weighted Alternative Least Squares) approach. However, the 
model-based methods are in general iterative methods, which 
usually converge after several iterations. 

A regular arrangement of station elements has the ad- 
vantage that it provides redundant baselines, i.e. baselines 
with the same physical length and orientation. Using redun- 
dant baseline informat i on for calibration was introduced by 
iNoordam and de Bruvnl (Il982l) . Its linearity, independence of a 
sky model, low computational cost and proven efficiency for 
precision calibration of WSRT (Westerbork Synthesis Radio 
Telescope) observations motivated us to apply the redundancy 
calibration to phased array stations. The redundancy calibration 
algorithm uses the data of all redundant baselines to obtain a 
convergent calibration solution in a single step. 

However, redundancy calibration in phased array systems re- 
quires additional considerations. This is essentially because of 
their different design concepts e.g. the closely located elements 
of a phased array station experience mutual coupling between 
elements which leads to non-identical beams of the station ele- 
ments and to correlated receiver noise. 

In this paper, we refine the standard data model presented 
in the phased array signal processing literature to account for 
baseline dependent corruptions in terms of the coupling effects. 
Using this refined data model, we briefly introduce the two cal- 
ibration methods i.e. model based and redundancy. This helps 
us to achieve a better understanding of the potential and lim- 
itations of both calibration methods. We wi ll also revise th e 
redundancy me thod formalism presented by dWieringal (ll99ll) : 
Liu et al. (2010)) to capture the nature of baseline dependent er- 
rors which affect the calibration accuracy. Some implementation 
issues will be raised and investigated using observed and sim- 
ulated data of LOFAR HBA stations. We will evaluate the re- 
dundancy calibration performance by comparing the variance of 
its results with the CRB (Cramer-Rao Bound) and the plots of 
residuals for the corrected data after redundancy calibration. We 
will also discuss limiting factors for its applicability. 

Although we have used HBA data to demonstrate the appli- 
cability and efficiency of the redundancy method, the analysis 
in this paper will be relevant to any phased array which is to be 
calibrated using redundant visibility information. 

Notation: we denote vectors in bold lowercase letters and 
matrices in bold uppercase letters. The matrix transpose and 
Hermitian transpose are denoted by (.) T and (.) H respectively. 
Operator diag(.) creates a diagonal matrix of a given vector. 
Operator vec(.) creates a vector. Operator G denotes element- 
wise matrix multiplication. Operator Z(.) returns the angle or 
phase of a complex number. Operator k(.) returns the condition 
number of a matrix. 



Map 1rom sky scan started at 14:10:20 UTC on 29 June 2009 
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Fig. 2: The sky imaged by HBA tiles at 14:10:20 UTC on 29 
June 2009. The Galactic plane appears in the north-west, the sun 
appears in the south-west. One can also see their corresponding 
grating response in the image. The image is presented in (I, m)- 
coordinates. I = cos{el)sin{az) and m = cos{el)cos(az), where el 
and az denote elevation and azimuth respectively. 



pling effects, a phased array of p elements has a signal vector, 
x(f) = [x\(t), X2(t), x p (t)] T which can be expressed as: 



x(f) = ro 



+ n(f) = TOAs(f) + n(f) 



(1) 



where s(f) is a q x 1 vector containing q mutually independent 
i.i.dQ Gaussian signals impinging on the array with the covari- 
ance of H s = diag(er s ) (size q x q), where cr s is a vecor of the 
source fluxes. They are also assumed to be narrow band, so we 
can define the q spatial signature vectors a* which include the 
phase delays due to the geometry and the directional response of 
the receiving elements (assumed to be identical). The vectors a* 
are called the array response vectors which are usually normal- 
ized. The receiver noise signals n,(f) are assumed to be mutually 
independent i.i.d. Gaussian signals in a p X 1 vector n(f) and 
uncorrected. Thus, the noise covariance matrix, E n = diag(<x„) 
(size p x p). The amplitudes and the phases of direction indepen- 
dent complex gains (g,'s) which have to be calibrated are y = 
[7i,72. —,7p] T an d </> = [e* , e-*, ...,e^ p ] T . Correspondingly 
r = diag(y) and O = diag(0). A = [ai.aa, ...,a 9 ] (size p x q) 
is a stack of the array response vectors. Before computing the 
coherency, x(f) is sampled with period T. The nth sample of the 
array signal vector x[n] is given by: 



x[n] = V x(t)6(t - nT) = x(nT) 



(2) 



2. Methods 

2.1. Data model for phased arrays 

The standard data model for phased array stations presented in 
the literature assumes that in the absence of RFI and any cou- 



N samples can be stacked in a matrix X = 
[x[l],x[2], ...,x[AH] (size p x N) which denotes the short 
term integrtion data set or snapshot. The array covariance matrix 
or the visibility matrix describing the correlation between all 



temporally independent and identically distributed. 
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sampled volatges can be estimated by R = XX 7 " /N whose 
expected value becomes: 

R = ra>Ai, s A H <3> H r H + E n (3) 

or, 

R = GAE S A H G H + £ n (4) 

Correlator errors can be represented as an additive term in the 
covariance matrix i.e. as a non-diagonal matrix. However in Eq. 
[3] and throughout this paper, they are disregarded assuming that 
the correlator is designed perfectly. 

Note that y,- could accommodate the overall amplitude gain 
of both the receiver system and the atmospheric disturbances and 
(pi the corresponding phase shift. In the case of station calibra- 
tion, we do not calibrate for direction dependent effects because 
we use snapshot data. In model-based methods, these are ab- 
sorbed in the known sky. In the case of redundancy calibration, 
these are absorbed in its fundamental assumption. We will elab- 
orate on this assumption in Sect. 12.31 

In phased arrays such as the HBA stations of LOFAR, or 
EMBRACE, the antenna elements are closely packed. This may 
cause mutual coupling between them, i.e. not all the power re- 
ceived by the elements is absorbed but some of the power is 
reradiated to the other elements. The reradiated power induces 
new currents in the other elements. Consequently, the radiation 
pattern of the elements changes. This leads to non-identicalness 
of the element beam patterns or different array response vec- 
tors despite their physical identicalness (as it is considered in 
A). This effect can be modeled as a direction dependent gain, Go 
(size pxq) being element-wise multiplied by the array response 
vectors: 

R = G(G O A)£ S (G () A) H M H G H + £ n (5) 

As it was studied bv lSvantessonl dl998l) . mutual coupling can 
directly be represented by M (size p x p) in the data model: 

R = GMLAX S A H M H G H + £ n (6) 

One may notice that a direct association of M and Go can not be 
expressed analytically. It requires a numerical evaluation which 
we present in Sect l3.ll For simplicity in order to continue our 
argument, we use the data model in Eq. [6] 

Mutual coupling may not only act on the signal, but also on 
the system noise. The LNA's (Low Noise Amplifiers) connected 
to the antennas in such an array, generate EM noise waves to- 
wards their outputs, but also send EM noise waves back into the 
antenna array. These waves are coupled into other receiver chan- 
nels, giving rise to a correlated noise contribution. This effect is 
known as noise coupling which contributes as R rec in a general 
and non-diagonal noise correlation matrix represented as: 

= H-sp + Rs/ty + R-ec (7) 

where R sp is the spillover noise correlation matrix which can 
usually be ignored as compared with R^ which is the sky noise 
contribution. Crosstalk in the back-end adds another baseline 
specific correlated noise term which we disregard here. Then, 
the general data model for the visibility matrix becomes: 

R = GMAE S A H M H G // + £; (8) 



Note that U a is not a diagonal matrix, unlike £ n . It can be 
shown that noise from each station element has still a normal 
distribution in view of large number of samples and according 
to the central limit theorem. Therefore, the off-diagonal elements 
of "L' n which will appear in Eq.[10]have a Wishart distribution. 

Although the noise model in Eq. [7] has been presented only 
in liter ature conce r ning P AFs (Phased Array Feeds) , for exam- 
ple by iJeffs et alj d2008l) and ilvashina etal] d2011l) it can be 
used a s a generic noise model for any antenna system (Maaskant 
d20 101) 1. It is clear, however, that depending on the antenna and 
LNA designs, the electrical characteristics of the array and the 
sparsity and density of the station elements layout, the strength 
of mutual coupling varies. Thus, the term R rec may be replaced 
by uncorrected receiver noise in this definition. 

2.2. Model-based calibration method 

The model based calibration problem has been for- 
mulated as a least squ ar es mi nimization problem by 
Wiin holds and van der Veen] d2009al) . Here, we rewrite it 
with the refined data model given in Eq|8] 

{g, <x n } = argmin g ^ ||GMAi: s A H M H G H + 2£ - R\\ 2 F (9) 

This estimates the noise and complex gain of each receiver 
element using the measured visibility, R and the modeled visibil- 
ity, GMA£ S A H M H G H + E s and A are assumed to be known. 
We can calculate them, if we specify the time of observation, the 
telescope geometry and known source parameters. In the pres- 
ence of the coupling effects (M and TS a ), estimation results are 
biased, unless an accurate beam model of each individual ele- 
ment is provided or the matrix M is known. In that case, one can 
handle the effect of the correlated noi se using a WALS approach 
(Wiin holds and van der V een ( 2009b])). 

Fig.|2]shows a sky map scanned by HBA tiles. One can see 
the Galactic plane due to the many short baselines and the sun as 
the dominant radio source. Given the beam models of the station 
elements, the WALS method treats the extended structure as cor- 
related noise and estimates the parameters of interest. However, 
the model based methods are in general iterative methods which 
usually converge after several iterations. 

2.3. Redundancy calibration method 

The basic assumption of redundancy calibration is that the re- 
dundant baselines sample the same Fourier component of the sky 
(AT, S A H ). This assumption holds if the array response vectors of 
the redundant baselines are the same i.e. the element beams are 
identical. We therefore begin with the most general data model 
given in Eq.[8]to understand the limitations of this method for a 
phased array station. 

To build up the system of equations for the redundancy cali- 
bration algorithm, we represent an off-diagonal element of R in 
Eq.[8]as: 

r,j = g ig ) [mae s a"m"]. . + z;. (io) 

Note that correlator offsets could contribute as additive cor- 
rupting factors in H n . . Since H s in Sect. 12. H is a diagonal matrix, 
one can expand Eq. [lOlto: 
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true visibility 



\Rij\ and if/jj are the amplitude and the phase of a complex visi- 
bility. Since we have to specify the absolute flux level, we set: 



Rij = g i g t j M ii M , jj Y u Ai q o- il A* 



(11) s 7i = o 



3=1 



I 

q=\ 



Yj Yj M ^ M *h» A ^' A M (r 'i 



Pl=\ p 2 = l 



additive term due to mutual coupling 



This shows that the presence of the mutual coupling pro- 
duces baseline dependent multiplicative and additive terms by 
influencing signals directly and noise indirectly (R rec in Eq. |7). 
This violates the fundamental assumption of redundancy i.e. we 
will not observe redundant visibilities on physically redundant 
baselines (see Fig. [4] top). 

Without loss of generality, MuMj. = 1 (or they can be ab- 
sorbed in the gains). To establish an analogy between Eq.fTTIand 
the redundancy method formalism in Wierin gal(ll99ll) . we define 
the following term: 



■ij = gig* x 

9=1 



P P 



Pl = l p 2 = l 

<.pi*i p 2 ?j 

and rewrite the Eq.fTTIas: 



+ 21 



a,ij 



(12) 



nobs _ ntrue * , 

R U ~ R U 8igj + e,j 



or 



R ™ gig , 



(13) 



(14) 



where R obs and R' rue are the observed and the theoretical redun- 
dant visibilities, respectively, g, and gj are the element complex 
gains. Wij can be defined as a baseline dependent error which 
affects the accuracy of the calibration results. 

We take the natural logarithm of both sides of the equation 
Eq.[14]to obtain: 



In R^ = InM + HYj\ + K9j ~ <M + In R'[" e + ln(w y ) 



(15) 



In the absence of mutual coupling, the first term of en drops 
and the second term drops the contribution of R rec . The only cor- 
related contributing terms come from R^y and R sp which are 
negligible depending on the SNR of the observation. Lets as- 
sume that we have such a case where we can ignore en. Then, 
we equate the amplitude and the phase values to decouple the 
system of equations for the phases and the amplitudes: 



ln\R°t s \= 7i + y } + W% 



Wij = <f>j ~ </>i + Wij 



(16) 



(17) 



(18) 



We also have to constrain the element phase. We can enforce this 
constraint by specifying that the average phase for all elements 
is zero: 



Zcpi = 



(19) 



Furthermore, there might also be an arbitrary linear phase slope 
over the array. This phase slope corresponds to a position shift 
of the field. This arises because a redundancy solution does not 
provide an absolute position. This can either be absorbed in the 
true visibilities or in the element phases. For a two dimensional 
array, we constrain x and y in the same manner: 



% t <l>iXi = 



Wi = 



(20) 



(21) 



where x, and yi are the (x,y) coordinates of the array elements. 
Eq.fTBT Eq. |2T| formulate the redundancy calibration method as 
two overdetermined systems of linear equations for phases and 
amplitudes which can be solved in a single step least square so- 
lution. For instance, the phase estimator can be symbolized as 
EphO - x P" bs and solved using the pseudo-inverse: 



6 



KM K^ obs 



(22) 



where = [<pi,(f> 2 , -, 4> P , WT"> ^2™"' •••> ^m"^ is the vector of P a " 
rameters to be estimated, m is the number of distinct redundant 
baselines, E p /, is the coefficient matrix and *P oil is the vector of 
the observed redundant phases and right sides of EqJT9l Eq]2T1 
Setting a phase reference for element phases, is done after the 
phase estimation. 

Note that we came to the solution in Eq. [22] because we ig- 
nored the baseline dependent noise. In its presence, our problem 
gets the form of E ph + = x ¥ obs (expanded in Eq. \ZQ. Then, 
the estimated parameters will deviate as follows: 



E ^e-e = \E T ph E ph \ l E T pl fi 



(23) 



where /? = vec(Z(ln(w, / ))). If the array mutual coupling is sig- 
nificant, the redundancy method will not be a reliable estimator. 
In the case of a weakly coupled array, the vector j8 still carries 
the correlated noise due to R sp and R^.. Assuming that these 
have Wishart distribution in Eq{12] their statistical distribution 
changes in vector j8. This has been taken into account for the 
results in Sect. 13.21 

Since the created systems of equations are highly sparse, 
they are computationally fast. Most importantly, they are inde- 
pendent of a sky model but their accu racy is aff e cted b y the SNR 
of the observed sky as discussed by iLiu etail (|2010). They re- 
quire identical beams of station elements and a regular arrange- 
ment of antennas to provide a sufficient number of redundant 
baselines. Then, one can solve for the element complex gains in 
a single step estimation. 
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3. Implementation of the redundancy calibration 
using real and simulated HBA data 

3. 1. Verification of the fundamental assumption of 
redundancy 

As mentioned earlier, redundancy translates to checking 
the identicalness of the element beams in a station. We 
have checked this for a 24-tile HBA station like the one 
shown in FigQ] using CAESAR (Co mputationally A d vance d 
and Efficient Simulator f or ARrays) dMaaskant et all (I2008I) : 
Maas kant and Yangl (l2006l) ). Numerical computation for the EM 
simulation was done using the CBFM (Characteristic Basis 
Function Method) which is the most numerically efficient 
and accurate method available for la rge sc attering problems 
dPrakash and Mittral (l2003l) : lYeo et al l (I2003I) ). Fig. [3] presents 
the simulation results when the beams are formed toward the lo- 
cal zenith. Each subplot presents the beam patterns of 24 tiles 
for a particular frequency in the = plane. One can see 
that due to a different mutual coupling environment, each tile 
has a slightly different beam pattern. The mutual coupling ef- 
fect is sufficiently small not to disturb the main beam but only 
the far sidelobes which are at least 12dB lower than the main 
lobe. However, it is not a favorable condition for redundancy 
calibration in general. Based on this simulation, we expect to 
observe non-redundant visibilities on redundant baselines when 
a strong source falls in the sidelobes. This has been confirmed 
by the real observation presented in Fig. |4] top. In this obser- 
vation, due to the absence of a strong source in the main beam, 
the non-identical sidelobes have the chance to play a significant 
role in disturbing the redundancy while in the observation shown 
in Fig.Hl bottom, one strong source in the main beam seems to 
be sufficient to dominate the influence of other possible strong 
sources which are observed through the non-identical sidelobes. 
There is of course a different contribution from the correlated 
background noise or sky noise (R^) in Fig. |U top and bottom 
panels. Fig.|4]presents two data sets obtained to verify the use of 
redundancy in a HBA station. The observations were made on 5 
September 2009. The tile beams in a HBA station (in the LOFAR 
core, known as CS302) were formed toward the local zenith. We 
let the sky drift over the FoV (Field of View) of the station. We 
captured the station visibilities approximately every 10 minutes 
(integration time of one second per frequency channel for 512 
frequency channels). The bandwidth of the frequency channels 



is 195 kHz and the frequency range is 100 - 200 MHz. On the 
right panel, the local sky viewed from CS302 is presentedn to 
show which sources have contributed to the visibilities shown 
on the left panel. 

Since, the mutual coupling environment changes at lower el- 
evations, the identicalness of the main beams is disturbed con- 
sequently up to 1 - 2dB. Real observations have shown less dis- 
turbance in the main beam at lower elevations (see Sect. l3.2l and 
Appendix. lAl. 

Based on the EM simulation by CAESAR (Fig. [3]) and 
real observations (Fig. |4}, we conclude that the best case sce- 
nario for redundancy calibration of an HBA station is to have 
a strong source in the main lobe when the beams are identical. 
However, the non-identicalness of the sidelobes introduces non- 
redundancy or systematic errors which can not be eliminated 
by any statistical method or longer integration time. Therefore, 
we must investigate whether the contribution of other strong 
sources, which may fall in the non-identical sidelobes, is sig- 
nificant. To quantify this systematic error on the visibility mea- 
sured on a given baseline, we present the following example. 
At 21:29:04 UTC on 24th November 2009, at an HBA station 
called RS208 located at (Ion = 6.9196°/a? = 52.6699°), four 
strong sources are in the FoV as shown in Fig. In Fig. the 
tile beams have been formed toward 3C461 or CasA. Thus, it 
is in the phase center. The sources 3C405 (CygA), 3C400 and 
3C144 (TauA) have fallen in the sidelobes. The tile beams and 
their standard deviation at 120MHz, at these sources are also 
depicted in Fig. The total complex visibility observed on a 
certain baseline, D,j is computed as: 



(25) 



where A(cr) = A(cr)/Ao is the normalized tile reception pattern 
at cr, Ao being the response at the beam center, /(er) is the source 
flux. Since there are four dominant strong sources in the FoV, 
we assume that the integral can be replaced by a summation. 
The complex plane in Fig. [6] right, shows summations of the 
visibility vectors as they were observed through hypothetically 
identical sidelobes, R foM / (solid line) as well as their summation 
when they are attenuated differently by the actual non-identical 
sidelobes, R totcll (dotted line). We compute the systematic errors 
in the phase and the amplitude of the visibility as follows: we 
compute the complex vector, R toto / = R3C461 + R3C144 + R3C405 + 



5 



P. Noorishad et al.: Redundancy Calibration of Phased Array Stations 



at MOO MHz 



atf=120 MHz 



atf=140 MHz 




20 40 60 80 
Theta (degrees) 

atf=160 MHz 



20 40 60 80 
Theta (degrees) 

atf=180 MHz 



20 40 60 
Theta (degrees) 

at f=200 MHz 




20 40 60 BO 

Theta (degrees) 



20 40 60 

Theta (degrees) 



20 40 60 ao 

Thela (degrees) 



Fig. 3: Simulated beam pattern of the tiles in a 24-tile HBA station using CAESAR. Each subplot presents the beam patterns of 24 
tiles in a particular frequency in <p — plane. The beams are formed toward the local zenith. The non-identicalness in the sidelobes 
is caused by the different mutual coupling environments of each tile. 



R3C400 using Eq. [25] For this, we choose the baseline type given 
in Fig. [8] right most panel (D) and frequency, / = 120MHz 
for = Til A, A is the wavelength. We set the averaged beam 
pattern of all 48 tiles as the identical beam for the term A((x) 
whose values are shown in Fig. [5] at different source locations, 
as the quantity beam level. We use A(&) - A((f)/A(<T3C46\)- 
The source fluxes, I(<f) are given in Table. Q] We also compute 

Kmi = R 3C46i + R 3ci44 + R 3C405 + R 3C40()- The computation set- 
tings are as before except that A(cr) is not identical this time. It 
is deviated from the averaged beam pattern by the standard devi- 
ation values shown in Fig. [5] at different source locations, as the 
quantity std(Beam). The bias introduced in the amplitude and 
the phase of R foMi as compared with R, OM ; are presented in the 
first row of Table, m 

We steer the beam toward each availabe strong source and 
repeat the calculation to predict the systematic error due to the 
other ones in the sidelobes with different standard deviations. 
The results are presented in Table. [1] We can conclude that track- 
ing CasA provides the best redundancy in the observed visibili- 
ties. 

The plot of the residuals for the corrected redundant visi- 
bilities shown in Fig. [8] confirms the results in Table. Q] Fig. [8] 
reveals non-redundancy in the measured visibilities (in terms of 
amplitudes and phases). In Table. [1] we predict the same quanti- 
ties which are caused by non-identical sidelobes. There is a few 
percent of discrepancy between the predicted residuals and their 
actual values in our observation. This can be explained by our 
simplifying assumption that we made for the number of strong 
sources in the FoV. However, Table. [T]can play an instrumental 
role for station calibration. 



3.2. Redundancy calibration performance 

To study the redundancy calibration performance, we tracked 
CasA on 24th November 2009 from 15:25:43 UTC until 
22:12:19 UTC. Running redundancy calibration on 48 captured 
data sets during this observation, gives very stable results for 
the receiver complex gains over time. This is an indication of 
system stability and a working calibration routine, which was 
running approximately every 10 minutes. We compare the vari- 
ance of the estimated complex gains over 48 runs of redundancy 
calibration with the CRB in Fig. [7] These quantities for the am- 
plitudes and the phases are presented in two separate plots as 
their estimators were decoupled in Sect. 12.31 The CRB (or the 
theoretical minimum variance) and the actual variance on the es- 
timated parameters over time are in good agreement. The small 
difference between them can be explained by not having exactly 
the same SNR from one observation to another during our sur- 
vey from 15:25:43 UTC until 22: 12: 19 UTC while the predicted 
S NR for this observation is SNR ^ 0.7 5, following the analysis 
bv lWiinholds and van Cappellenl d201 ll) . This value was used to 
compute the presented CRB in Fig. [7] 

We also studied the residuals for the corrected redundant vis- 
ibilities. Fig.[8]shows an example of the residuals in amplitudes 
and phases of the corrected visibilities on a distinct type of re- 
dundant baseline. The snapshot is captured at 21:29:04 UTC, 
when CasA is at high elevation. The first row shows the results 
after redundancy calibration as compared to the results after the 
model based calibration in the second row. The type of redun- 
dant baseline is depicted on the station configuration in the right 
most panel. The integration time is 1 second per frequency chan- 
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Fig. 4: An observation with a HBA station on September 5, 2009 at 13:13:48 UTC in which redundancy calibration fails (top) and 
an observation with the same sation at 23:02:01 UTC in which redundancy calibration is successful (bottom). The left panels show 
the measured visibilities for three distinct type of redundant baselines. The right panels show the local sky model at the time of 
observation with strongest sources on the sky superimposed on a contour plot of the element beam pattern at 170 MHz. 

Table 1: A-team radio sources in the FoV of RS208, on 24th of November 2009 at 21:29:04 UTC and the predicted levels of 
the systematic errors in the amplitudes and the phases of the redundant visibility due to the contribution of other strong sources 
through non-identical sidelobes. The baselines given in Fig. [8] right most panel and frequency, / = 120MHz were chosen for the 
computations. 



Calibrator source 


source flux [Jy] 


normalized error in amplitudes 


error in phases [rod] 


3C461 (CasA) 


8609 


0.07 


0.06 


3C405 (CygA) 


8100 


0.20 


0.09 


3C400 


540 


0.93 


0.01 


3C144(TauA) 


1420 


0.84 


0.28 



nel. More plots of the residuals from the same observation are 
presented in AppendixlAl Fig. IA. H and Fig. IA.2I By comparing 
them, one may notice that the baseline length does not make a 
significant difference in the residuals. The residuals are in the 
order of 2 - 5% in both phases and amplitudes. Because a strong 
source like CasA in the field of view dominates the effect of pos- 
sible correlated noise, this can be explained by not having 100% 
identical main beams or the corruption due to the nonidentical 



sidelobes (see Table. [T). However these results are satisfactory. 
Moreover, Fig.|7]showed, the estimated complex gains which are 
our parameters of interest in a station calibration, are very sta- 
ble. This is due to the constraints on them in the formalism of 
the redundancy calibration. 

Using the data captured at 15:25:43 UTC, when CasA is at 
low elevation, reveals slightly larger residuals in the order of 
5 - 10% (see AppendixjA] Fig. IA.31 Fig. IA.5l ). This is due to a 
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Fig. 5: A-team radio sources in the FoV of an HBA station called RS208 on 24th November 2009 at 21:29:04 UTC. Tile beams are 
streered toward 3C461 to provide the most redundant visibilities. The values of the averaged reception pattern of 48 tiles, A(<x) and 
their standard deviation, std(A{cr)) at 120MHz have also been indicated at other sources locations, cr. 



different mutual coupling environment which leads to less iden- 
tical beams either in the main lobe or the sidelobes. 

The residuals also show that the two calibration methods per- 
form almost equally well, although redundancy method behaves 
more consistently. Model based calibration has slightly larger 
residuals on short baselines, e.g. up to ~ - 1% in phases and 
~ 0-2% in amplitudes whereas, it shows similar residuals as re- 
dundancy calibration on long baselines. Extended structures e.g. 
the galactic plane or north polar spur are captured on short base- 
lines. Modeling them is computationally expensive and some- 
how impractical. Therefore, in the model based method, one has 
to discard the visibilities measured on short baselines to simplify 
the measured sky for a corresponding simple sky model while 
the redundancy method is sky model independent. Experiment 
has shown that discarding the equations of short baselines in re- 
dundancy calibration routine, will not improve its results signif- 
icantly. 

Moreover, model based methods are sensitive to RFI sources, 
as their presence confuses the sky model. Since redundancy cal- 
ibration is independent of the sky model, it is less sensitive to 
RFI. One may have noticed that more frequency channels had to 
be flagged for model based calibration. 

As discussed by iLiu et al.l (HoiO), redundancy calibration 
quality depends on the SNR of the observed sky although it is 
independent of a sky model. ILiu et al .1 (120101) showed that the 



estimated parameters are affected differently in the presence of 
baseline dependent noise assuming it is Gaussian noise. The vec- 
tor fS in Eq.l23lreveals a non-Gaussian baseline dependent noise. 
We have taken this into consideration for the following results. 
After adding different levels of a non- i.i.d. Gaussian noise to the 
output vector of the array (Eq. Q3, we simulated the station vis- 
ibility assuming that the station elements have identical beams. 
We used the station configuration of RS208 (see Fig. [8] right 
most panel). The complex gains were estimated over 100 runs 
of Monte-Carlo simulation for different SNR per visibility. The 
variance of the error in the estimated complex gains (over sim- 
ulation runs with different random noise) versus different SNR 
per visibility is shown in Fig. [9] These plots indicate the vari- 
ance of the estimated complex gains which are our parameters 
of interest in a station calibration, are the same for the differ- 
ent tiles. This shows the stability of the algorithm. However the 
analysis like this is very array configuration dependent, as the 
array configuration determines the coeficient matrices by which 
our estimators are defined. 

3.3. Limitations of the redundancy calibration 

The applicability of redundancy calibration is limited by the fol- 
lowing factors: 

1 . The station configuration 
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Fig. 6: Total complex visibility on a certaion baseline due to the contribution of all strong sources shown in FigJ5] The ones which 
contribute through non-identical sidelobes introduce non -redundancy or systematic errors to the total redundant visibility. This 
has been demostrated in a complex plane on the right. The solid line shows the visibility vectors as they were observed through 
hypothetically identical sidelobes. The dotted line shows the visibility vectors when they are attenuated differently by the actual 
non-identical sidelobes. 
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Fig. 7: Theoretical minimum variance of the estimated amplitudes and phases of the complex gains or CRB compared with their 
actual variance. These quantities for the amplitudes and the phases are presented in two separate plots as their estimators were 
decoupled in Sect. 12.31 



2. Mutual coupling between the station elements which leads 
to: 

- non-identical element beams 

- presence of baseline dependent noise 

3. SNR of the observed sky 



The station configuration can influence the suitability of re- 
dundancy in different ways. The redundancy calibration method 
requires a regular arrangement of station elements. In a station 
with p elements one needs a sufficient number of distinct types 
of redundant baseline to have a system of equations in which 



9 



P. Noorishad et al.: Redundancy Calibration of Phased Array Stations 



0,2 




ItlO l?0 1*5 ISC ISO ?oo 

freq [MHz] 



T3 




UM l?0 140 1jtS HO XIO 

freq [MHz] 












illi 






IfWf 







& LI.043 



'WWR... 




freq [MHz] 



freq [MHz] 



pn ir. 1<1 



Fig. 8: Plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant base- 
lines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibration. 
The data is taken from the observation done on 24th November 2009 at 2 1 :29:04 UTC, when CasA is at high elevation. The station 
configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility are 
also, with the same color code depicted on the station layout. 




all station element gains are involved. The more redundant the 
station's baselines are, the less information in the measured vis- 
ibilities will be missed in the computation. This is not a concern 
in the HBA stations, as they are highly redundant. We define the 



ratio of the number of the measured visibilities that are used in 
the computations to the total number of the measured visibili- 
ties as /. The station configuration also determines the coeficient 
matrices i.e. E ; ,/, or E amp i for phase and amplitude estimators 
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respectively. The condition number is a relative error magnifi- 
cation factor i.e. errors in the right-hand side of a linear system 
of equations can cause errors k(E p /,) times as large in the solu- 
tion. Fig. [10] shows three different configurations for the HBA 
stations in the LOFAR system. Given a certain configuration, af- 
ter recognizing the redundant baseline in it, one can calculate the 
aforementioned quantities. These are given in Table. [2] This ta- 
ble shows that RS208 is the most reliable configuration, from the 
redundancy calibration point of view. However, these quantities 
can be figures of merit when we design for redundancy within 
the stations and in arrangement of the stations within the whole 
array. 

Configuration and element spacing in a station are usually 
decided based on obtaining low sidelobes of the station beam. 
Minimizing the mutual coupling between the station elements 
should be taken into consideration as other important figures of 
merit in a station layout for the sake of redundancy applicabil- 
ity. This consideration serves a twofold goal, obtaining identical 
element beams at least in the main beam and minimizing the re- 
ceiver correlated noise (R )ec ). Non-identicalness of the element 
beams introduces systematic errors which can not be eliminated 
by any statistical method or longer integration time. In the pres- 
ence of the baseline dependent noise, the estimated parameters 
will deviate as given by Eq. [23] Quantifying this deviation for 
a given array, requires a good understanding of the noise terms 
in Eq. [7] especially R rec . The correla ted noise terms are hard to 
model analytically, Maaskant] d2010l) . One may use a powerful 
and generic tool called CAESAR to compute them numerically 
using its EM- (ElectroMagnetic) and MW- (Micro Wave) simula- 
tors. This is an ongoing work for HBA and EMBRACE stations. 
Primary results from the MW-simulation for an HBA station has 
shown a negligible contribution of R rec . However, Fig. [9] indi- 
cates that whatever that effect on calibration accuracy is, it does 
not make redundancy algorithm unstable. 

4. Discussion 

In this paper, we studied the applicability of redundancy calibra- 
tion to phased array stations for the first time. Its performance 
was demonstrated using data of a new telescope, LOFAR. It re- 
quired new considerations which were not part of the original 
design of conventional arrays e.g. WSRT and VLA. We took 
them into account by refining the data model. Reformulating 
the redundancy calibration formalism using the new data model, 
helped us to understand its potential, limitations and the effects 
of non-Gaussian baseline dependent noise on the calibration ac- 
curacy. 

Obtaining identical beams for station elements and minimiz- 
ing the correlated noise, have to be considered as figures of merit 
for the configuration and spacings between the elements within 
a station in SKA (Square Kilometre Array) pathfinders and the 
SKA itself. This is definitely necessary, regardless of the cal- 
ibration method to be used. For the redundancy method, it is 
fundamental, while having identical element beams saves com- 
putational capacity and time for model based methods. 

In the plots of residuals and Table. [TJ we demonstrated how 
non-identicalness of the beams leads to systematic errors which 
can not be eliminated by any statistical methods or longer inte- 
gration time. Since, identical element beams are fundamental es- 
pecially for the redundancy calibration, similar EM-simulations 
to what we presented here, are highly recommended for any ar- 
ray which is to be calibrated using redundancy. Moreover, re- 
dundancy turns out to be an extremely suitable diagnostic tool 
for recognizing and monitoring failing station elements. This is 



of concern in arrays with large number of elements such LOFAR 
and SKA. It is currently being used for this purpose in the 
LOFAR system. 

Wide field science wit h new radio telescopes such as SKA 
(van Ardenne et al . (2009)), demands immense signal process- 
ing and computational capacity. This is mainly due to the large 
number of elements, variable beams and ionosphere over a large 
field of view. The redundancy calibration results encourage us to 
recommend to seriously consider redundancy in the SKA config- 
uration, both at station and at whole array level. In large interfer- 
ometers such as SKA, where the baselines are non-coplaner, one 
has to think of redundancy in UVW-space, however, not only in 
UV-space. Redundancy can be applied for gain calibration of the 
whole array in specific regimes where the stations observe the 
sky, through the same ionospheric patch. It also saves computa- 
tional capacity and gives more accurate estimate of the telescope 
gains as compared to the model based gain calibration. This is 
of course a major step forward to achieve radio images of higher 
dynamic range. 

5. Conclusion 

We studied for the first time the applicability and limitations of 
redundancy calibration in phased array stations using both real 
and simulated data of a radio telescope using the new aperture ar- 
ray technology, LOFAR. The results clearly show that the addi- 
tional contraints from redundant baselines do improve the qual- 
ity of the calibration and in addition provide a powerful tool for 
system diagnostics at different levels of the telescope phased ar- 
ray hierarchy (both intra station and inter station). The merit of 
redundancy in the station and full array layout and the additional 
advantage for diagnostic are well demonstrated by this study. 
We therefore strongly recommend to design for redundancy in 
both the station layout and the array configuration of future aper- 
ture arrays, in particular the SKA, where the required dynamic 
range will be an order of magnitude beyond any existing array. 
Because redundancy gives a better handle on characterizing the 
state of the system, it provides in addition a model independent 
diagnostic tool for subsystems such as a station. 
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Fig. 10: Three different station configuration for HBA stations. 
Table 2: Station configuration and the condition of redundancy calibration 



Station Conf. 


no. of distinct redundant baseline 


size(Ep/,) 


size(E am ,,,) 




K(Eampl) 


I 


RS208 


84 


1127 x 132 


1125 x 132 


52.72 


9.74 


99.65% 


CS302 


113 


1123 x 161 


1121 x 161 


6.74el6 


9.99 


99.29% 


CS004 


72 


547 x 120 


545 x 120 


8.42el6 


1.73el5 


48.22% 
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Appendix A: More plots of residuals for corrected 
redundant visibilities 

Some more plots of the residuals for the corrected redundant vis- 
ibilities of our observational campaign on 24th November 2009 
are presented here. Fig. lA.ll and Fig. lA.2l are similar results of the 
same data set (captured at 21:29:04 UTC, when CasA is at high 
elevation) presented in Fig|8]but for different baselines, shown 
on the right most panel. 

Fig. IA.31 Fig. IA. 51 show similar results using the data cap- 
tured at 15:25:43 UTC, when CasA is at low elevation. They 
reveal slightly larger residuals in the order of 5 - 10% as com- 
pared with the results in Fig. [8] Fig. IA. 1 1 and Fig. IA.2I This is 
due to a different mutual coupling environment which leads to 
less identical beams either in the main lobe or the sidelobes. As 
it was discussed in Sect. 13.21 we still see smaller residuals af- 
ter redundancy calibration as compared with the model based 
method. 
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Fig. A.l: Plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant 
baselines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibra- 
tion. The data is taken from an observation done on 24th November 2009 at 21:29:04 UTC, when CasA is at high elevation. The 
station configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility 
are also, with the same color code depicted on the station layout. 
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Fig. A. 2: Plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant 
baselines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibra- 
tion. The data is taken from the observation done on 24th November 2009 at 21:29:04 UTC, when CasA is at high elevation. The 
station configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility 
are also, with the same color code depicted on the station layout. 
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Fig. A. 3: The plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant 
baselines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibra- 
tion. The data is taken from the observation done on 24th November 2009 at 15:25:43 UTC, when CasA is at low elevation. The 
station configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility 
are also, with the same color code depicted on the station layout. 
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Fig. A. 4: The plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant 
baselines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibra- 
tion. The data is taken from the observation done on 24th November 2009 at 15:25:43 UTC, when CasA is at low elevation. The 
station configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility 
are also, with the same color code depicted on the station layout. 
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Fig. A. 5: The plots of the residuals for corrected redundant visibilities in terms of phases and amplitudes on a given set of redundant 
baselines. The first row shows the results after redundancy calibration. The second row shows the result afters model based calibra- 
tion. The data is taken from the observation done on 24th November 2009 at 15:25:43 UTC, when CasA is at low elevation. The 
station configuration of RS208 is shown in the right most panel. The corresponding redundant baselines to each redundant visibility 
are also, with the same color code depicted on the station layout. 



17 



